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ABSTRACT 


A technique  of  statistically  adjusting  dynamical  fore- 
casts of  tropical  cyclone  motion  was  tested.  All  tests  were 
performed  with  operationally-analyzed  data  from  the  U.S.  Navy 
Fleet  Numerical  Weather  Central  (FNWC).  Three  sets  of  regres 
slon  equations  were  developed  to  modify  forecasts  of  typhoon 
tracks.  The  first  set  of  equations  was  based  only  on  forward 
Integration  of  the  FNWC  Tropical  Cytlone  Model  (TCM)  for  28 
cases  In  1975-76.  An  independent  sample  of  cases  from  1977- 
78  indicated  that  the  first  equation  set  was  based  on  too 
small  a sample  size,  especially  considering  the  anomalous 
nature  of  the  1975-76  storm  tracks.  A second  equation  set 
based  only  on  forward  integration  of  the  TCM  was  derived  from 
61  storm  track  forecasts  from  1975-78.  Results  from  the  ex- 
periments with  these  equations  indicate  that  systematic  data 
and  model  errors  can  be  used  to  statistically  adjust  forecast 
storm  tracks.  The  second  equation  set  based  on  forward  inte- 
gration showed  improvement  over  the  unmodified  model  predic- 
tions at  all  forecast  times.  A third  equation  set  based  on 
forward  and  backward  integration  of  the  TCM  explained  the 
greatest  amount  of  variance  of  all  the  equation  sets.  In  a 
dependent  test  of  these  equations  using  31  of  the  1977-78 
cases,  the  U.S.  Navy  7th  Fleet  error  goal  of  100  and  150 
nautical  miles  at  48  and  72  hours  was  nearly  met. 
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I.  INTRODUCTION 


The  tropical  storm  Is  one  of  nature's  most  destructive 
phenomena,  making  typhoon  prediction  a major  concern  world- 
wide. The  need  to  forecast  storm  conditions  has  produced  a 
variety  of  prediction  methods.  Subjective  methods,  includ- 
ing persistence,  have  generally  produced  fairly  good  short- 


range  forecasts.  An  especially  difficult  problem  for  any 
subjective  or  objective  prognostic  scheme  is  the  occurrence 
of  storm  recurvature.  This  problem  was  significant  in  the 
western  North  Pacific  Ocean  during  the  1975  typhoon  season, 
and  to  a lesser  extent  during  the  1976  typhoon  season.  An 
unusual  number  of  storms  tracked  northward,  while  others  re- 
curved. As  a result,  the  forecast  error  statistics  were 
higher  than  normal  (Annual  Typhoon  Report,  1975).  The 
United  States  Seventh  Fleet  Commander,  noting  the  necessity 
to  improve  the  forecast  errors,  has  levied  a requirement  for 
the  Joint  Typhoon  Warning  Center  (JTWC)  to  achieve  maximum 
forecast  errors  of  50,  100  and  150  nautical  miles  for  24, 

48  and  72  hours,  respectively.  In  comparison,  the  1977  JTWC 
average  position  forecast  errors  for  tropical  cyclones  at  24-, 
48-,  and  72-hours  were  140,  266  and  390  nautical  miles,  re- 
spectively (Annual  Typhoon  Report,  1977).  The  long-term 
error  trend,  as  indicated  by  the  five-year  mean  error,  has 
been  increasing  since  1972.  In  1977  the  operational  TCM 
produced  mean  vector  errors  of  138,  262  and  450  nautical 
miles  for  all  tropical  cyclones  (Annual  Typhoon  Report,  1977). 


Although  progress  has  been  made  In  the  numerical  simula- 
tion of  tropical  storm  characteristics  such  as  Intensity  and 
spiral  bands  (Anthes  et  al  , 1974),  the  major  emphasis  has 
focused  on  storm  motion  (Hovermale  et  al  , 1976  ; Ley,  1975). 
Multi-level,  nested-grid  models  are  being  developed  In  an 
attempt  to  Improve  forecast  positions  over  those  of  analog, 
persistence  and  statistical  techniques.  In  spite  of  the 
sophistication  of  some  numerical  models,  undesirable  conse- 
quences may  result  from  Initialization  with  poor  or  limited 
data.  One  test  of  the  dynamical  models  Is  the  accuracy  of 
the  Initial  storm  track.  If  the  short-term  forecasts  are 
Inaccurate,  a dynamic  model  cannot  be  expected  to  produce 
accurate,  extended  forecasts.  Some  dynamic  models  (Hovermale 
et  al  , 1976  ; Ley,  1975;  Hodur  and  Burk,  1978)  predict  tracks 
with  systematic  bias  relative  to  the  actual  track  and  com- 
monly predict  motion  which  is  too  slow.  A likely  source 
of  error  seems  to  be  in  specification  of  the  Initial  data 
fields,  although  inadequate  resolution  in  the  numerical 
model  could  also  explain  the  biases. 

Since  the  resources  required  to  increase  high  quality 
data  or  to  run  sophisticated  numerical  models  are  expensive 
and  time  consuming,  alternative  methods  for  improving  model 
forecasts  should  be  considered.  Statl sti cal -dynami cal 
schemes  for  predicting  tropical  storm  motion  (Renard  et  al, 
1972)  rely  heavily  on  current  storm  motion.  If  position 
reports  of  tropical  storms  are  not  timely,  then  the  statisti- 
cal schemes  may  be  Impaired.  An  approach  which  circumvents 
many  of  these  problems  Is  to  use  the  model  forecast  positions 


themselves  to  statistically  modify  the  predicted  track  posi- 
tions. It  Is  the  primary  objective  of  this  thesis  to  de- 
velop and  evaluate  statistical  regression  equations  for 
adjusting  dynamically  predicted  storm  tracks  from  the  Tropi- 
cal Cyclone  Model  (TCM)  used  by  the  U.S.  Navy  Fleet  Numerical 
Weather  Central  (FNWC). 

The  basic  model  used  for  these  experiments  is  the  primi- 
tive-equation, three-layer,  tropical  cyclone  model  developed 
by  Elsberry  and  Harrison  (1971)  and  Harrison  (1973).  Although 
this  model  is  capable  of  triply-nested  operation  (Harrison, 
1973),  results  from  Ley  and  Elsberry  (1976)  show  that  the 
coarse  and  nested  grids  produced  nearly  identical  results  In 
a selected  case  study  based  on  hand-analyzed  data.  In  1975 
FNWC  adapted  a coarse-mesh  (2°),  three-layer,  dry  version  of 
this  model  for  tropical  cyclone  prediction.  Preliminary 
results  with  operational  data  were  presented  by  Hlnsman 
(1977)  for  1975  and  1976  data.  For  the  1977  typhoon  season, 
the  model  was  modified  to  include  a biasing  technique  sug- 
gested by  Shewchuk  and  Elsberry  (1978).  The  use  of  the 
forecast  stream  function  field  to  determine  the  cyclone 
forecast  positions  reduced  errors  in  the  relative  vorticlty 
tracking  used  prior  to  this  modification  (Shewchuk,  1977). 

The  current  operational  version  of  the  model  has  boundary 
conditions  which  are  Insulated,  free-sllp  walls  on  the 
north  and  south  and  cyclic  on  the  east  and  west.  Since  it 
was  suspected  that  these  boundary  conditions  could  adversely 
affect  the  forecast  storm  track,  Hodur  and  Burk  (1978)  In- 
corporated one-way  interactive  boundaries.  Substantial 
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Improvement  In  storm  track  prediction  was  made  compared  to 


the  previous  channel  version  of  the  model, 
boundary  version  of  the  TCM  was  used  In  the 
described  herein. 


The  open 
experiments 


II.  BACKWARD  INTEGRATION 


A.  THE  MODEL 

. A detailed  description  of  the  Tropical  Cyclone  Model 

' l 

'i  (TCM)  Is  presented  In  Appendix  A.  The  TCM  Is  a three-layer, 

primitive-equation  model  In  pressure  coordinates.  It  Is 
an  open  boundary  model  with  one-way  Interactive  boundary 
conditions  on  the  north  and  south  and  cyclic  conditions  on 
the  east  and  west  boundaries. 

■ 

B.  MODEL  MODIFICATIONS 

All  forward  Integration  of  the  model  was  carried  out 
using  the  version  of  the  TCM  described  by  Hodur  and  Burk 
(1978).  The  backward  Integration  was  carried  out  using  a 
negative  time  step  of  -600  seconds  unless  the  northern 
boundary  location  was  greater  than  40°N.  In  the  latter 
case,  the  negative  time  step  was  reduced  to  -450  seconds 
after  the  storm  moved  a sufficient  distance  to  the  north. 

In  the  forward  Integration,  heat  was  added  to  the  storm 
center  as  defined  by  a minimum  wind  at  1000  mb.  The  pur- 
pose of  the  heating  function  Is  to  counteract  the  dispersion 
of  the  vortex  due  to  the  finite  differencing  (Ley  and 
Elsberry,  1976).  If  heat  had  been  added  during  the  back- 
ward Integration,  this  would  have  contributed  to  a better 
definition  of  the  storm  center.  However,  this  presents  a 
physically  unrealistic  situation  for  a typhoon  moving  back- 
ward  In  time.  For  this  reason,  the  heating  function  was 
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set  equal  to  zero  for  all  backward  Integrations.  Negative 
heating  was  rejected  because  It  led  to  premature  dispersion 
of  the  storm,  which  became  difficult  to  track. 


III.  APPROACH  TO  TRACK  MODIFICATION 


Two  types  of  regression  equations  were  tested  here.  One 
type  will  have  predictors  based  only  on  forward  Integration 
of  the  TCM , while  the  other  set  of  predictors  will  be  based 
on  both  forward  and  backward  Integration  of  the  TCM.  The 
open  boundary  version  of  the  TCM  (Hodur  and  Burk,  1977;  Hodur 
and  Burk,  1978)  was  run  with  46  cases  (9  storms)  from  1975 
and  1976.  All  cases  were  based  on  operationally-analyzed 
data  obtained  from  FNWC.  The  regression  equations  developed 
from  this  sample  were  then  tested  against  Independent  cases 
In  1977  and  1978.  A second  setof  regression  equations  based 
only  on  forward  Integration  of  the  TCM  are  derived  using  the 
combined  data  sets  of  1975-78. 

The  final  set  of  regression  equations  are  based  on  forward 
Integration  of  the  model  to  72  hours,  as  well  as  backward 
Integration  for  36  hours.  The  backward  integration  should 
reveal  the  effects  of  systematic  model  and  data  errors. 
Fundamental  assumptions  are  that  the  model -rel ated  errors 
tend  to  be  systematic,  and  adjustment  for  data  errors  is 
possible  where,  In  the  absence  of  observations,  the  Initial 
analysis  reverts  to  the  east-west  flow  appropriate  to  the 
climatology.  Backward  Integration  will  increase  the  number 
of  predictors  available  to  explain  the  variance  between  TCM 
forecasts  and  best  track  positions.  This  should  lead  to 
Improved  regression  equations  with  a higher  explained  variance. 
The  predictors  using  backward-integrated  positions  can  be 
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used  to  adjust  the  forward-integrated  positions  using  the 
same  Initial  data. 

Both  the  forward-  and  backward-predicted  tracks  were 
compared  to  the  best  tracks  at  their  corresponding  times. 
Next;  regression  equations  to  adjust  the  predicted  track  to 
the  actual  track  using  stepwise  regression  were  generated 
with  the  Statistical  Package  for  the  Social  Sciences  (SPSS). 
The  regression  equations  were  developed  to  make  corrections 
every  12  hours  up  to  72  hours.  If  a TCM  run  was  Incomplete, 
alternate  sets  of  equations  were  derived  according  to  the 
length  of  the  TCM  run.  The  options  considered  were  36-,  48-, 
60- , and  72-hour  TCM  runs,  with  and  without  backward  Integra- 
tion. 

If  successful,  the  advantage  of  this  approach  would  be 
that  use  of  simple  regression  equations  would  require  much 
less  computer  time  than  more  sophisticated  dynamic  models. 
Thus,  It  may  be  possible  to  produce  tracks  that  are  more 
accurate  than  -the  ordinary  open  boundary  TCM  without  a 
large  Increase  In  computer  resources.  An  operational  advan- 
tage of  this  approach  Is  that  no  warning  positions  are 
required  for  computation  of  any  predictors  In  the  regression 
equations.  The  method  and  results  of  this  approach  to  storm 
track  prediction  are  presented  In  the  following  sections. 


I ' 
I i 


IV.  DEVELOPMENT  OF  THE  REGRESSION  EQUATIONS 

A.  THE  PREDICTORS  ANO  PREDICTANOS 

Using  TCM  forecast  runs  versus  best  tracks  (Annual  Typhoon 
Reports,  1975-77),  12  predictands  were  derived  by  computing 
the  east-west  and  north-south  differences  between  positions  at 
corresponding  times  (see  Fig.  1).  Storm  positions  at  each 
forecast  time  are  adjusted  by  two  regression  equations,  one 
for  the  east-west  direction  and  the  other  for  the  north-south 
direction.  Thus,  for  a 72-hour  TCM  forecast  run,  a total  of 
12  regression  equations  would  be  used  to  modify  storm  tracks 
in  12-hour  increments  from  12  to  72  hours. 

Predictors  used  in  these  equations  were  model -predi cted 
displacement  and  velocity,  broken  into  components  along  the 
east-west  and  north-south  directions.  The  Julian  day,  lati- 
tude and  longitude  of  the  initial  position  of  each  TCM  fore- 
cast run  are  also  included  as  predictors  in  the  regression 
equations.  A schematic  illustration  of  the  intervals  over 
which  the  predictors  were  calculated  for  the  forward  integra- 
tion runs  is  shown  in  Fig.  2.  A complete  list  of  predictands 
and  predictors,  with  the  times  for  which  they  were  computed, 
appears  in  Table  I. 

B.  METHOO  OF  EQUATION  DERIVATION 

The  regression  equations  were  derived  using  the  Statisti- 
cal Package  for  the  Social  Sciences  (file,  et  al,  1975). 

Cases  with  missing  values  for  predictors  or  predictands  were 
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Figure  1. 


Depiction  of  the  model  errors  (i.e.,  the  predic- 
tands),  which  are  the  difference  between  the  best 


track  and  forecast  positions,  which  are  shown 


above  as 


etc . 


OA  » 


Figure  2.  Depiction  of  intervals  over  which  displacements 
and  average  speeds  were  computed  using  only  for- 
ward integration  of  the  open  boundary  TCM.  Un- 
modified TCM  forecast  positions  (0)  are  shown  at 
12  hour  Intervals. 


TABLE  I 


L 


Predl ctors/predi ctands  used  to  develop  regression  equations 
for  typhoon  track  modification  based  only  on  forward  Inte- 
gration of  the  open  boundary  TCM. 


Predictands:  AX,  AY 

Times  at  which  predictands  are  computed: 

12,  13,  26,  38,  60,  72  hrs  , 

Predictors:  AX,  AY,  u,  v 

Time  intervals  over  which  each  predictor  was  computed 


00-12, 

12-24, 

24-36,  36-48, 

48-60, 

60-72 

00-24, 

12-36, 

24-48,  36-60, 

48-72, 

00-36 

0k 

00 

1 

o 

o 

00-60, 

00-72  hrs 

Initial  Position  Predictors:  Julian  Day,  Latitude, 

Longi tude 


Ti mes  : 

- Initialization  time  of  a given  TCM  forecast  run 
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automatically  eliminated  from  all  calculations.  Such  cases 
arise  because  the  tracking  routine  is  not  always  able  to 
follow  the  storm  center  throughout  the  72-hour  interval.  In 
the  forward  integration  tests  of  Hodur  and  Burk  (1978)  with 
the  open-boundary  conditions,  only  28  of  the  46  cases  extended 
to  72  hours. 

To  avoid  the  problem  of  uncomputable  predictors  when 
the  duration  of  TCM  runs  was  less  than  72  hours,  it  was 
decided  to  derive  alternate  sets  of  equations  based  on  the 
duration  of  the  TCM  forecast  run.  The  forecast  lengths  con- 
sidered were  36,  48,  60  and  72  hours.  For  example,  if  the 
TCM  produced  only  a 48-hour  storm  track  forecast,  no  regres- 
sion adjusted  track  would  extend  behond  48  hours,  and  only 
predictors  in  the  range  00-48  hours  were  considered  when 
deriving  the  regression  equations  to  be  applied  to  a 48-hour 
TCM  run. 

All  predictors  listed  in  Table  I were  considered  for  in- 
clusion in  each  regression  equation.  Selection  of  predictors 
was  stopped  when  the  next  variable  in  the  stepwise  regression 
explained  less  than  1%  of  the  variance. 

C.  FIRST  EQUATION  SET  BASED  ON  FORWARD  INTEGRATION 

The  1975-76  cases,  as  well  as  those  in  1977  and  1978  to 
be  discussed  later,  are  based  on  cases  with  relatively  well- 
developed  storms  which  are  more  readily  modeled  by  the  coarse- 
mesh  TCM.  It  should  be  noted  that  the  operational  version  of 
the  TCM  is  used  only  for  storms  exceeding  50  knots.  Inclusion 
of  weaker  tropical  storms  would  likely  increase  the  variance 


between  the  forecast  and  best  track  positions,  thus  making 
It  more  difficult  to  derive  stable  regression  equations.  A 
sample  equation  for  the  regression  adjustment  along  the  x- 
axis  at  72  hours  is  shown  below: 

DXER 72=*-537. 015 1+39. 2088 (XXLAT) +68. 1605 (VX1 224) 

- 50. 1251 (VX0024)-43.8154(VY0048) 

- 1.4956(0X6072)+  0 . 6 102 ( JULDAY ) 

Velocity  was  the  most  frequently  selected  predictor.  This 
result  is  not  too  surprising  if  one  recalls  that  a common 
fault  of  dynamic  models  is  prediction  of  motion  which  is  too 
slow  (Hovermale  et  al , 1976;  Ley  and  Elsberry,  1976).  The 
version  of  the  TCM  (Hodur  and  Burk,  1978)  used  in  these  ex- 
periments is  no  exception.  Predictors  with  their  associated 
time  intervals  beginning,  ending  or  overlapping  the  valid 
time  of  a given  regression  equation  were  often  selected.  In 
general,  this  suggests  a sort  of  "statistical  extrapolation" 
in  which  past,  present  and  future  forecast  motions  are  used 
to  correct  the  forecast  track. 

The  average  explained  variance  of  the  regression  equations 
appears  in  Table  II.  This  parameter  measures  the  strength 
of  the  linear  relationship  between  the  multi-linear  regres- 
sion equation  value  and  the  observed  value.  The  amount  of 
explained  variance  in  Table  II  generally  decreases  as  the 
TCM  runs  become  shorter  in  duration  and  the  number  of  avail- 
able predictors  is  thereby  reduced.  In  general,  longer  TCM 
forecast  runs  permitted  use  of  regression  equations  which 
corrected  for  more  of  the  variance  from  the  best  track. 
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TABLE  II 


Average  explained  variance  of  the  regression  equations 


Forward  Integration  (1975-76  Cases) 


72  Hr  TCM  Run 
60  Hr  TCM  Run 
48  Hr  TCM  Run 
36  Hr  TCM  Run 


81.6 * 
77.6 * 
75.4* 
68.4* 


Forward  Integration  (1975-78  Cases) 


72  Hr  TCM  Run 
60  Hr  TCM  Run 
48  Hr  TCM  Run 
36  Hr  TCM  Run 


51.5* 

46.5* 

47.6* 

43.2* 


Forward  and  Backward  Integration  (1977-78  Cases) 


72  Hr  TCM  Run 
60  Hr  TCM  Run 
48  Hr  TCM  Run 
36  Hr  TCM  Run 


86.2* 

84.3* 

83.8* 

82.8* 


Shortcomings  (see  Fig.  3)  of  the  one-way  (OW)  interac- 
tive boundary  version  of  the  TCM  are  that  the  forecast 
track  typically  runs  to  the  left  of  the  storm  path  and  the 
velocities  are  generally  too  slow  (Hodur  and  Kurk,  1978). 
These  systematic  errors  of  the  TCM  seem  to  contribute 
significantly  to  the  large  amount  of  explained  variance  in 

the  regress  ion-produced  typhoon  positions  (see  Table  II). 

/ 

It  should  be  recalled  that  the  1975  typhoon  season, 
and  to  a lesser  extent  the  1976  season,  experienced  a high 
frequency  of  storm  recurvature  as  well  as  a large  number  of 
storms  which  tracked  northward.  This  bias  was  present  in 
the  data  sample  used  to  derive  the  regression  equations  and 
had  considerable  impact  on  the  selection  of  predictors  and 
the  computation  of  constants  and  coefficients. 

D.  SECOND  EQUATION  SET  BASED  ON  FORWARD  INTEGRATION 

To  Increase  the  number  of  typhoon  cases  used  to  derive 
the  regressio-n  equations,  TCM  forecast  runs  from  1977  and 
1978  were  included  in  the  data  set.  These  TCM  forecasts 
were  originally  used  to  make  an  independent  test  of  the  1 975- 
76  regression  equations.  However,  the  majority  of  the 
typhoons  in  1975  and  many  in  1976  had  significant  recurvature 
The  inclusion  of  the  1977-78  storm  tracks  should  make  the 
sample  more  representative  of  the  various  tracks  found  in 
the  western  North  Pacific  Ocean  area. 

Because  some  TCM  runs  did  not  extend  to  72  hours,  only 
61  of  90  cases  were  available  for  derivation  of  the  regres- 
sion equations.  As  with  the  first  set  of  equations  based 
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on  the  forward  Integration  of  the  TCM,  alternate  sets  of 
equations  were  derived  based  on  72-hour,  60-hour,  48-hour 
and  36-hour  TCM  runs.  When  less  than  1%  of  the  variance  was 
explained  by  a variable,  selection  of  predictors  was  halted. 
The  complete  set  of  regression  equations  Is  listed  In 
Appendix  B. 

The  expanded  data  sample  (61  cases)  contained  a greater 
variety  of  storm  tracks  and  reduced  the  recurvature  bias  of 
the  1975-76  typhoon  cases.  Storm  tracks  In  the  1975-78 
data  set  used  In  these  experiments  are  characterized  by  four 
general  categories:  westward,  northwestward,  northward, 
and  recurving  paths.  With  this  greater  variety  of  storm 
paths,  the  amount  of  explained  variance  of  the  regression 
equations  (see  Table  II)  dropped  commensurately . For  example, 
the  regression  equations  for  a 72-hour  TCM  run  incurred 
approximately  a 30  percent  decrease  In  explained  variance. 
With  the  exception  of  a small  variation  at  48  hours  in  Table 
II,  the  general  trend  was  again  a reduction  of  explained 
variance  as  the  duration  of  the  TCM  forecast  runs  decreased 
to  36  hours.  As  discussed  previously,  this  characteristic 
Is  attributed  to  having  fewer  predictors  available  to  explain 
the  variance  when  the  TCM  forecast  is  of  shorter  duration. 
These  results  reaffirm  that  this  method  of  adjusting  TCM 
forecasts  Is  at  Its  best  when  the  TCM  forecast  extends  to 
72  hours,  rather  than  a shorter  forecast  Interval. 

Velocity  was  again  the  most  frequently  selected  type  of 
predictor.  Indicating  an  attempt  to  compensate  for  the  over- 
all slowness  of  the  TCM.  However,  there  was  a better 
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balance  of  velocity  and  displacement  predictors  In  the  1975- 
78  equations  as  compared  to  the  1975- 76  equations. 

E.  EQUATION  SET  BASED  ON  FORWARD  AND  BACKWARD  INTEGRATIONS 

The  third  set  of  regression  equations  was  based  on  storm 
positions  derived  from  both  backward  and  forward  Integration 
of  the  TCM . Additional  velocity  and  displacement  predictors 
based  on  backwards  Integration  of-  the  TCM  as  Indicated  In 
Fig.  4 and  Table  III  were  added  to  the  data  set.  These  pre- 
dictors are  computed  In  the  same  manner  as  those  based  on 
forward  Integration.  All  predictors  listed  In  Table  I and 
Table  III  were  considered  for  Inclusion  In  each  regression 
equation.  Regression  equations  were  again  developed  for  TCM 
forecast  runs  of  36,  48,  60  and  72  hours  duration. 

The  regression  equations  which  included  backward-inte- 
grated positions  had  the  highest  average  explained  variance 
(see  Table  II).  As  TCM  forecast  runs  become  shorter  In 
duration,  the  amount  of  explained  variance  in  the  equations 
decreases  due  to  fewer  predictors  being  available  (see 
Table  II).  The  complete  set  of  regression  equations  Is 
1 i sted  i n Appendix  B . 

The  regression  equations  took  advantage  of  the  systematl 
errors  inherent  In  the  initial  fields  and  In  the  numerical 
model.  The  most  favored  predictors  from  the  backward  Inte- 
gration occurred  In  the  Interval  from  -12  to  00  hours.  For 
the  12  equations  used  when  the  TCM  was  Integrated  forward 
to  72  hours  and  backward  to  -36  hours,  predictors  In  the 
Interval  -12  to  00  hours  appeared  in  9 equations.  The 


Figure  4.  Depiction  of  intervals  over  which  displacements 
and  average  speeds  were  computed  using  both  for 
ward  and  backward  integration  of  the  open  bound 
ary  TCM.  Unmodified  TCM  forecast  positions  (0) 
are  shown  at  12-hour  intervals. 


TABLE  III 

Adcltlonal  predictors  used  to  develop  regression  equations 
for  typhoon  track  modification  based  on  the  backward  Inte- 
gration of  the  open  boundary  TCM. 


Predictors:  AX,  AY,  u,  v 

Times  * 00-M1 2 , M12-M24,  M24-M36,  12-M12, 

00-M24 , M12-M36 , 00-M36,  24-M24  hrs 


whe  re 


00  * initial  time 

M24  = minus  24  hours  from  initial  time 
12  a plus  12  hours  from  Initial  time 
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velocity  along  the  y-axls  from  -12  to  00  hours  appeared  In  7 
of  the  equations , and  In  each  case  It  explained  the  most 
variance  of  any  predictor  In  the  equation.  Approximately 
35  percent  of  all  the  predictors  In  these  12  regression 
equations  were  computed  from  backward-integrated  storm 
positions.  Velocity  was  the  most  frequently  selected  pre- 
dictor from  all  forecast  Intervals.  As  will  be  Indicated 
In  later  examples,  the  regression  equations  appeared  to  be 
compensating  for  the  generally  slow  motion  of  the  TCM.  It 
should  be  noted  here  that  these  regression  equations  and 
the  sample  cases  which  follow  were  based  on  a limited  data 
set  of  31  storm  cases  from  1977  and  1978. 
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V.  RESULTS 


A.  FIRST  EQUATION  SET  BASED  ON  FORWARD  INTEGRATION 

The  primary  purpose  of  the  tests  with  the  1975-76  cases 
was  to  make  a preliminary  evaluation  of  the  regression 
equations  which  were  derived  using  this  data  set.  Because 
this  data  set  was  used  to  derive  the  regression  equations, 
these  experiments  are  referred  to  as  a dependent  test. 

A sample  of  adjusted  storm  positions  based  on  the  first 
equation  set  is  shown  in  Figs.  5-7.  Results  with  Typhoon 
June  shown  In  Fig.  5 and  Typhoon  Marie  depicted  in  Fig.  6 
were  very  encouraging.  Typhoon  Marie  is  the  same  storm 
which  appears  in  Fig.  4.  Storm  recurvature  was  accurately 
depicted,  and  storm  velocity  was  markedly  improved  in  these 
two  storms.  Note,  however,  that  the  velocities  in  the 
regression  tracks  are  generally  greater  than  the  best  track 
velocities.  This  characteristic  of  the  regression  equations 
appears  to  be  an  attempt  to  compensate  for  the  slowness  of 
the  TCM.  In  the  majority  of  cases  in  this  sample,  the 
storms  were  excessively  accelerated  due  to  this  feature  of 
the  equations.  Since  the  regression  equations  were  based 
primarily  on  predictors  derived  from  TCM  forecasts,  it  is 
clear  that  the  goodness  of  the  regression-adjusted  storm 
track  is  dependent  on  the  quality  of  the  TC.M  forecast  itself. 
The  adjusted  track  for  Typhoon  Rita  (August  1975)  shown  in 
Fig.  7 illustrates  how  a poor  TCM  forecast  will  lead  to  an 
extremely  radical  regression  adjustment.  This  example 
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Figure  6.  Typhoon  Marie  (10  Apr  76,  00  GMT),  otherwise  the 
same  as  Fig.  5. 


Figure  7.  Typhoon  Rita  (21  Aug  75,  12  GMT),  otherwise  the 
same  as  Fig.  5. 


Indicates  how  small  TCM  velocities,  as  well  as  radical 
changes  In  the  forecast  storm  track,  will  result  In  a 
regression-adjusted  track  with  extremely  high  velocities 
and  unrealistic  variations  In  the  storm  path.  This  behavior 
of  the  regression  equations  seems  to  indicate  that  obviously 
erroneous  regression  positions  can  be  used  as  a basis  to 
reject  TCM  forecast  positions  as  well. 

The  mean  forecast  errors  of  the  1975-76  dependent  test 
sample  are  listed  in  Table  IV.  The  sample  improved  on  the 
forecast  errors  of  both  JTWC  and  the  unmodified  TCM  fore- 
casts for  1975-76.  At  72  hours  with  a sample  size  of  28,  the 
regression  equation  errors  were  also  less  than  the  U.S. 

Navy  7th  Fleet  goal  of  150  nm.  However,  it  is  not  expected 
that  this  pattern  would  be  repeated  in  an  independent  sample 
of  typhoon  cases. 

B.  INDEPENDENT  TEST  OF  FIRST  EQUATION  SET 

An  independent  sample  of  44  cases  from  the  1977-78  typhoon 
seasons  was  then  used  to  test  the  regression  equations  de- 
rived from  the  1975-76  cases  discussed  in  the  previous 
section.  As  indicated  by  the  statistics  in  Table  V,  the 
results  of  the  independent  test  were  generally  very  poor. 

In  each  category  the  modified  TCM  tracks  were  worse  than  the 
unmodified  tracks.  The  large  errors  are  attributed  to  the 
unstable  nature  of  the  regression  equations  from  the  small 
and  very  homogeneous  sample  of  anomalous  storm  tracks  used  to 
derive  the  regression  equations.  This  led  to  excessive 
velocities  and  erratic  tracks  in  the  regression-adjusted 
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storm  positions  for  1977-78.  When  the  TCM  forecasts  had 
large  deviations  from  the  best  track,  the  regression  positions 
Incurred  commensurately  larger  errors.  It  was  clear  that 
the  number  of  TCM  forecast  cases  had  to  be  Increased  as  much 
as  possible  to  Increase  the  stability  of  the  regression 
equations . 

The  regression  equations  used  In  these  tests  did  not 
Improve  storm  positions  by  taking  advantage  of  the  systematic 
errors  of  the  TCM,  which  is  often  too  slow  and  to  the  left 
of  the  best  track.  Typhoon  Lola  is  depicted  in  Fig.  8 and 
indicates  slight  corrective  shifting  of  the  regression  posi- 
tions to  the  right  of  the  TCM  forecast  track.  However,  over- 
compensation is  again  evident  in  the  velocity  components. 

More  typical  of  the  independent  test  results  is  the  adjusted 
track  for  Typhoon  Lucy  shown  in  Fig.  9.  Velocity  errors  were 
large  and  the  track  was  erratic.  In  many  cases  in  the  inde- 
pendent test,  an  erratic  track,  such  as  the  regression- 
adjusted  path  -for  Typhoon  Lucy,  occurred  when  the  unmodified 
TCM  forecast  track  also  incurred  large  errors.  The  single 
redeeming  feature  of  these  1975-76  equations  may  be  that 
highly  erratic  regression  positions  could  be  indicative  of  a 
poor  TCM  forecast.  This  information  may  be  useful  in  leading 
the  typhoon  forecaster  to  reject  both  the  modified  and  the 
unmodified  TCM  guidance  in  making  his  decision. 

C.  SECOND  EQUATION  SET  BASED  ON  FORWARD  INTEGRATION 

The  second  set  of  regression  equations  was  obtained  from 
the  combined  storm  cases  in  the  1975-76  and  1977-78  samples. 
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Figure  8.  Typhoon 
same  as 
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Because  this  exhausts  the  sample,  these  equations  can  only 
be  evaluated  as  a dependent  test. 

This  second  set  of  equations,  with  a larger  sample  of 
61,  was  more  stable  overall  than  the  first  set  of  regression 
equations.  While  the  regression-adjusted  positions  still 
increase  the  velocity  excessively,  the  apparent  result  is 
that  this  overcompensation  is  decreased  in  the  1975-78  re- 
gression equations  versus  those  for  1975-76. 

The  five  typhoon  cases  shown  previously  in  Fig.  5 through 
Fig.  9 appear  again  in  Fig.  10  through  Fig.  14,  but  this  time 
illustrating  the  regression-adjusted  positions  of  the  1975- 
78  equations.  The  track  for  Typhoon  June  shown  in  Fig.  10 
was  not  adjusted  as  well  in  this  test,  but  it  serves  to 
illustrate  that  the  regression  adjustments  cannot  always  com- 
pensate for  the  tendency  of  the  TCM  to  not  predict  recurvature. 
Presumably  if  the  TCM  track  had  been  more  northwesterly,  the 
adjustment  would  have  been  toward  more  recurvature.  The 
track  for  Typhoon  Marie  depicted  in  Fig.  11  still  indicates 
a reasonable  adjustment  for  recurvature,  but  has  incurred  a 
marked  decrease  in  velocity  adjustment.  Although  the  adjusted 
speed  of  movement  for  Typhoon  Rita  shown  in  Fig.  12  is  slower 
than  that  on  Fig.  7,  the  track  is  not  significantly  changed 
from  the  previous  result.  However,  this  case  again  suggests 
that  a radical  and  obviously  erroneous  regression  adjustment 
may  serve  as  a basis  to  reject  the  TCM  forecast  also.  Typhoon 
Lola  (September  1978),  illustrated  in  Fig.  13,  has  a much 
improved  regression  adjustment  of  the  storm  track  compared 
with  Fig.  8.  The  velocity  corrections  in  this  case  are  much 
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Figure  10.  Typhoon  June  (18  Nov  75,  00  GMT),  otherwise  the 
same  as  Fig.  5. 
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Figure  13.  Typhoon  Lola  (29  Sep  78,  00  GMT),  othe 
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REGRESSION  ADJUSTMENT 


closer  to  the  actual  velocities  which  appear  in  the  best 
track  positions.  The  adjusted  track  for  Typhoon  Lucy  is 
shown  in  Fig.  11  and  is  to  be  compared  with  Fig.  9.  The 
regression-adjusted  positions  are  still  poor,  but  not 
nearly  as  erratic.  Note  that  the  regression-modified  posi- 
tions do  suggest  some  recurvature,  although  inaccurately, 
and  that  the  decreased  velocities  are  much  more  realistic 
with  the  1975-78  regression  equations. 

Generally,  the  results  of  these  experiments  suggest 
that  a further  increase  in  the  number  of  TCM  forecasts  in 
the  data  sample  would  lead  to  a further  improvement  in  the 
regression-modified  storm  positions.  Overall,  the  errors 
incurred  by  the  regression-adjusted  tracks  were  less  than 
those  of  the  unmodified  TCM  storm  path  (see  Table  VI).  The 
improvements  are  especially  noteworthy  at  48  hours  and  72 
hours . 

In  forecasting  typhoons  which  tracked  westward,  the  un- 
modified TCM  forecast  positions  usually  proved  to  be  very 
accurate,  with  the  regression-modi fi ed  path  seldom  improving 
on  the  TCM.  One  case  (Typhoon  Rita,  October  1978)  in  which 
the  regression  equations  improved  the  forecast  is  shown  in 
Fig.  15.  The  improvement  is  the  result  of  an  Increase  in 
storm  velocity  by  the  regression  equations.  Several  other 
examples  of  enhanced  track  predictions  are  shown  In  Fig.  16 
through  Fig.  18.  These  examples  show  some  improvement  in 
direction,  but  the  Improvement  In  speed  of  movement  is 
especially  rewarding. 
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Fig. 


D.  EQUATION  SET  BASED  ON  FORWARD  AND  BACKWARD  INTEGRATION 
These  experiments  with  the  1977-78  typhoon  cases  were 
performed  to  test  the  regression  equations  derived  from  the 
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1977-78  data  set  which  included  backward-integrated  posi- 
tions to  -36  hours.  As  the  complete  set  of  backward  tracks 
was  used  in  the  derivation  of  the  regression  equations, 
these  experiments  refer  only  to  the  dependent  sample. 

The  results  of  the  dependent  test  with  forward-  and 
backward-integrated  positions  were  the  best  of  all  experi- 
ments conducted  during  this  research.  This  equation  set  in- 
curred the  lowest  mean  errors  overall  (see  Fig.  19).  The 
regression-modified  storm  position  errors  are  very  close 
to  meeting  the  U.S.  Navy  7th  Fleet  error  goal  (see  Table  VII). 

Tracks  for  Typhoons  Lola,  Rita,  and  Babe  shown  previously 
in  Fig.  13,  15,  and  17,  respecti vely,  appear  again  in  Fig.  20 
through  Fig.  22  with  adjusted  positions  based  on  forward- 
backward  integration.  Typhoon  Lola  (September  1978),  shown 
in  Fig.  20,  is' more  erratic  than  the  same  case  (Fig.  13) 
using  regression  equations  based  only  on  forward  integration. 
However,  the  regression  positions  that  are  based  on  forward- 
backward  integration  are  fluctuating  on  either  side  of  the 
best  track.  The  variation  about  the  best  track  is  attributed 
primarily  to  the  small  size  of  the  data  set  (31  cases)  used 
to  derive  the  equations.  This  "saw-tooth"  variation  about 
the  best  track  is  characteristic  of  the  behavior  storms  which 
tracked  westward  in  this  experiment.  As  in  previous  tests, 
it  proved  difficult  to  Improve  on  the  unmodified  TCM  forecast 
positions  when  typhoons  were  tracking  approximately  westward. 
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FOREC  ST  TIME  (HRS) 


Figure  19 


Mean  track  error  (nm)  for  the  1975-76  dependent 
test  (0),  the  1977-78  independent  test  (a),  the 
1 975-78  dependent  test  (•)  , and  the  1977-78 
dependent  test  (o) . The  1977-78  dependent  test 
(o)  was  the  only  experiment  which  included  back 
ward  integration. 


TABLE  VII 
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The  adjusted  track  of  Typhoon  Rita  in  Fig.  21  exhibits  a 
similar  behavior  In  another  westward  tracking  storm,  but 
does  very  well  In  the  longitudinal  positioning  of  the  storm. 
Adjustment  of  the  track  of  Typhoon  Babe,  shown  In  Fig.  22, 
has  Improved  the  velocity  forecast  as  well  as  has  indicated 
some  recurvature  which  the  TCM  missed  (see  also  Fig.  17). 

Typhoon  Phyllis  is  Illustrated  in  Fig.  23.  It  provides 
an  excellent  forecast  of  recurvature  not  achieved  by  the 
unmodified  TCM  forecast  or  by  the  regression  equations  based 
only  on  forward  integration  of  the  TCM.  The  velocities  in- 
dicate only  a slight  tendency  to  overcompensate  for  the 
slower  TCM  forecast.  Recurvature  of  Typhoon  Gilda,  depicted 
in  Fig.  24,  is  correctly  predicted,  but  the  velocities  are 
excessive.  Adjustment  of  the  track  of  Typhoon  Wendy  as  in 
Fig.  25  showed  slight  improvement  over  the  unmodified  TCM 
positions.  The  interesting  feature  in  the  regression-adjusted 
path  is  the  indication  that  the  storm  would  reverse  its 
direction.  Within  12  hours  of  the  indicated  time,  Typhoon 
Wendy  did  change  direction  in  this  manner,  but  with  somewhat 
smaller  displacements. 
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V I . CONCLUSION S 


A three-layer,  primitive-equation  model  (Hodur  and  Burk, 
1978)  with  one-way  interactive  boundaries  is  being  tested 
at  the  Naval  Environmental  Prediction  Research  Facility 
(NEPRF).  The  objective  of  this  research  was  to  generate 
statistical  regression  equations  to  adjust  the  TCM- predi cted 
tracks  towards  the  best  tracks.  This  approach  is  based  on 
the  assumption  that  it  is  possible  to  adjust  for  systematic 
model  and  data  errors.  Development  consisted  of  deriving 
three  sets  of  regression  equations,  with  two  sets  based  only 
on  forward  integration  of  the  TCM,  while  the  last  set  con- 
tained predictors  based  on  both  forward  and  backward  inte- 
gration. All  TCM  runs  were  based  on  operationally  analyzed 
data  from  FNWC.  Ti me- dep enden t boundary  conditions  provided 
to  the  TCM  were  derived  from  analyzed  (rather  than  predicted) 
fields,  as  was  the  case  in  Hodur  and  Burk  (1978). 

The  most  notable  improvements  occured  with  the  regression 
equations  containing  predictors  based  on  backward  integration 
of  the  TCM.  The  equations  with  predictors  having  both  for- 
ward- and  backward-integrated  positions  explained  the  greatest 
amount  of  variance  of  any  set  of  regression  equations.  Ad- 
justments to  the  TCM  tracks  at  12-hour  intervals  resulted 
In  predictions  in  the  1977-78  dependent  test  that  nearly 
met  the  7th  Fleet  forecast  error  goals.  The  selection  of  the 
velocity  along  the  y-axis  in  the  interval  -12  to  00  hours  as 
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the  predictor  explaining  the  most  variance  suggests  strongly 
that  systematic  errors  at  or  near  initialization  time  were 
used  advantageously  to  adjust  TCM  track  forecasts.  It  is 
noteworthy  to  observe  that  TCM  forecasts  with  all  data  samples 
used  in  these  experiments  incurred  systematic  errors  (see 
Fig.  26). 

It  should  be  noted  that  the  weak  link  in  the  regression 
adjustment  scheme  is  the  small  sample  size  for  derivation  of 
the  regression  equations.  With  the  exception  of  the  1375-76 
equations,  the  regression  coefficients  have  not  been  tested 
with  an  independent  sample  of  storms.  Application  of  the 
statistical  equations  in  the  future  may  not  produce  comparable 
results  due  to  the  relatively  small  number  of  cases  used  to 
derive  the  equations.  Likewise  some  years  have  persistently 
anomalous  storm  tracks,  as  in  the  first  sample  (1975-76) 
treated  here.  It  is  expected  that  some  increase  in  the 
stability  of  the  regression  equations  could  be  achieved  by 
increasing  the^  number  of  TCM  forecast  cases  used  to  derive 
the  equations.  This  is  especially  true  for  the  backward 
integration  set. 

In  all  of  the  experiments,  the  unmodified  TCM  forecasts 
were  difficult  to  beat  in  the  case  of  westward  propagating 
storms.  However,  occasional  improvement  in  velocity  predic- 

i 

tion  was  noted.  The  typhoons  used  in  these  tests  had  storm 
tracks  in  several  general  categories,  such  as  westward,  south- 
east-to-northwestward,  northward,  or  recurving.  Since  the 
regression  equations  were  derived  from  storm  cases  involving 
various  combinations  of  these  track  types,  it  was  felt  that 
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these  widely  varying  tracks  and  subsequent  errors  contri- 
buted to  failures  of  the  regression  equations  to  make  proper 
track  adjustments  In  some  cases.  A possible  solution  might 
be  to  vastly  increase  the  TCM  case  sample  size  and  then 
subdivide  the  cases  Into  “regression  analogs"  according  to 
the  direction  of  typhoon  propagation.  The  next  step  would 
be  to  derive  a set  of  regression  equations  based  on  each 
subset  of  storm  tracks.  It  could  tre  expected  that  each  sub- 
set of  equations  would  have  inherent  adjustment  biases 
towards  the  type  of  track  from  which  they  were  derived.  Any 
or  all  of  these  subsets  of  equations  could  be  applied  to  each 
unmodified  TCM  forecast. 

On  the  basis  of  the  sample  examined  here,  it  appears  that 
adjustment  for  recurvature  was  handled  most  effectively  by 
the  regression  equations  which  included  backward-i ntegrated 
tracks.  It  is  suggested  that  regression  equations  based  on 
an  expanded  sample  of  TCM  cases  including  backward  integration 
should  be  derived  to  produce  better  adjustments  and  statis- 
tical equations  with  greater  stability. 
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APPENDIX  A:  THE  TROPICAL  CYCLONE  MOOEL 


A.  THE  FORECAST  MODEL 

The  coarse-mesh  version  of  the  primitive  equation  model 
developed  by  Elsberry  and  Harrison  (1971),  Harrison  (1973), 
Ley  and  Elsberry  (1976),  and  Hodur  and  Burk  (1978)  was  used 
as  the  basis  for  these  experiments.  The  model  Is  a coarse 
mesh  (2°),  three-layer,  dry  model  with  one-way  interactive 
boundaries  on  the  north  and  south  and  cyclic  east-west 
boundaries.  The  equations  used  in  the  model  are: 
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L(S)  represents  the  flux  divergence  of  any  scalar  quantity  S. 
The  meteorological  symbols  used  above  can  be  found  In  the 
"List  of  Symbol s" . 

A sufficient  condition  for  the  linear  computational  sta- 
bility of  the  solution  for  two-dimensional  equations  govern- 
ing simple  wave  motion  is: 
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where 

C = phase  velocity  of  the  fastest  gravity  wave 
At  = time  increment 
Ax  = horizontal  grid  increment 
Computational  stability  requires  a maximum  time  steo  o 4*0 
seconds  for  this  model.  Ley  (1975)  achieved  an  increase  to 
800  seconds  by  time  averaging  the  pressure  gradient  term  of 
the  momentum  equations.  Shewchuk  (1977)  used  a 600  second 
time  step  for  testing  his  1975  cases  as  the  relocatable  grid 
was  extended  northward  where  Ax  is  reduced.  In  the  TCM  ver- 
sion (Hodur  and  Burk,  1978)  used  for  these  experiments,  the 
time  step  was  600  seconds  unless  the  northern  boundary  exceeded 
40°N,  in  which  case  the  time  step  was  450  seconds. 

The  initial  step  was  forward  in  time,  while  the  leapfrog 
time  differencing  scheme  was  used  in  all  subsequent  iterations. 
Friction  was  neglected,  and  the  consequent  storm  motion  was 
primarily  the  result  of  advective  processes  and  heating.  A 
Bessel  interpolator  was  used  to  locate  the  minimum  wind  at 
1000  MB.  Then  latent  heating  was  simulated  by  adding  heat  to 
a horizontal  7x7  grid  centered  on  the  minimum  wind,  which 
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was  not  necessarily  on  a grid  point.  Weighting  of  the  heat- 
ing function  In  the  vertical  was  0.3,  1.0,  and  0.3  at  850, 

550,  and  250  MB  respectively.  In  the  horizontal,  the  effect 
of  heating  Is  smoothed  out  by  a Cressman  weighting  function 
(Haltlner,  1971)  and  results  in  a less  erratic  storm  track. 

B.  THE  GRID 

The  forecasts  are  carried  out. on  a uniform,  coarse-mesh 

✓ 

C2°)  Mercator  grid  true  at  22.5°N.  The  horizontal  grid  Inter- 
val was  205.8  km.  The  domain  consisted  of  32  points  east- 
west  and  24  points  north-south.  The  grid  was  oriented  so 
that  each  storm  was  initially  located  southeast  of  the  center 
of  the  grid.  The  vertical  distribution  of  variables  is  shown 
in  Fig.  A-l.  Although  the  variables  are  staggered  in  the 
vertical,  staggering  on  the  horizontal  grid  was  not  used  in 
these  experiments.  The  850  mb  winds  were  used  to  compute  the 
advective  terms  at  1000  mb. 

C.  BOUNDARY  CONDITIONS 

Boundary  conditions  on  the  walls  were  after  Hodur  and 
Burk  (1978).  In  the  one-way  interactive  (QW)  model,  wind  com- 
ponents normal  to  the  boundaries  are  adjusted  so  there  is  no 
net  divergence  from  the  forecast  domain.  By  Integrating 

H'vn  <A-8> 

around  the  boundary,  the  corresponding  ip  values  are  obtained 
(Hawkins  and  Rosenthal,  1965).  Distance  along  the  boundary 
is  represented  by  s,  while  Vn  is  the  velocity  component  normal 


Fig.  A-l.  Vertical  distribution  of  dependent  variables 
and  pressure  levels  for  the  three-dimensional 
model  (after  Harrison,  1973). 
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to  the  boundary  with  the  positive  direction  being  Inward. 

Date  from  outside  the  forecast  grid  must  continuously 
be  Incorporated  into  the  boundaries.  This  is  accomplished 
as  described  by  Perkey  and  Kreltzberg  (1976),  and  following 
their  notation,  the  prediction  equation  for  a variable  X 
becomes 


xnd,j) 


X p ( 1 , j )+W( 1 , j ) 


3 x_ 
m 

Tr 


At+[1-W(i  , j ) ] 


1 . J 


3Xis 

TE 


At 


1,j 


( A — 9 ) 

where  the  subscripts  n and  p indicate  the  new  and  the  previous 
value  of  X.  The  subscript  m refers  to  the  model  tendency, 
while  Is  is  the  large-scale  tendency.  The  weighting  function 
used  in  this  model  was  somewhat  different  from  that  given  by 
Perkey  and  Kreitzberg  (1976). 

The  following  weighting  function  W(i,j)  was  constructed 
so  that  a minimum  amount  of  noise  was  produced  near  the 
boundaries 

0.0  on  the  boundaries 
0.05  one  grid  row  in  from  boundaries 
0.25  two  grid  rows  in  from  boundaries 
0.45  three  grid  rows  In  from  boundaries 
0.65  four  grid  rows  In  from  boundaries 
0.85  five  grid  rows  in  from  boundaries 
1.00  on  all  other  Interior  points 
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The  forecast  fields  produced  must  be  filtered  near  the 
boundaries.  The  filter  used  In  this  model  is 

FJ-  (l-lFj  ♦ Ct>(FJ+1  ♦Fj.,}  (A"  10) 

where  Fj  is  the  filtered  data  at  point  j,  F ^ is  the  unfiltered 
data  at  point  j,  and  a represents  the  smoothing  parameter. 

The  filter  was  applied  every  40  minutes  to  the  six  rows  and 
columns  nearest  the  boundaries  with  a = 0.5.  This  value  of 
a yields  a response  function  of  0 for  a 2Ax  wave.  This  filter 
is  also  applied  over  the  entire  grid  to  all  the  prognostic 
variables  every  three  hours.  In  this  manner  the  forecasts 
include  large-scale  tendencies  from  outside  the  forecast 
domain. 

The  vertical  velocity  at  the  upper  boundary  is  equal  to 
zero  and  Is  calculated  at  levels  5,  3,  and  1,  Fig.  A-l, 
through  downward  integration  of  Eq.  (A-5). 

In  an  operational  mode,  the  time  dependent  boundary  condi- 
tions must  be  specified  by  a global  forecast  model.  In  these 
experiments,  forecast  fields  were  not  available.  Therefore, 
using  a " perfect-prog"  approach,  analyzed  fields  taken  every 
12  hours  were  used  to  specify  the  boundary  values. 

D.  INITIALIZATION 

Analyzed  fields  obtained  from  FNWC  are  used  as  input  data 
for  the  model.  Initialization  of  the  model  is  accomplished 
by  calculating  non-di  vergent  winds  from  the  stream  function  ip . 
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The  relative  vortlclty,  the  forcing  function  for  the  stream 
function,  was  obtained  from  the  observed  u and  v components 

5 * m2i7x(FTJ  " ly^1  (A-ll) 

The  normal  component  of  the  wind  was  adjusted  so  there  Is  no 
net  Inflow  or  outflow  In  the  domain.  The  stream  function  was 
found  from  the  expression 


( A- 12 ) 


A direct  solver  (see  Faulkner  and  Rosmond,  1976)  rather  than 
successive  over-relaxation  is  used  to  efficiently  solve  the 
Poisson  equations  In  the  initialization  process.  The  non- 
divergent  wind  components  were  then  calculated  through 


u 


( A — 1 3 ) 


An  appropriate  balance  between  the  mass  and  motion  fields 
Is  achieved  through  partial  differentiation  of  Eq.  ( A - 1 ) with 
respect  to  x and  Eq.  (A-2)  with  respect  to  y.  Addition  of 
these  equations  and  assuming  that  the  time  rate  of  change  of 
divergence  and  gradients  of  the  map  scale  factor  can  be  neg- 
lected leads  to  a Poisson  equation 


V2<p 


, ~ « 3v . 3u  . 

' ff(I7CL(%)]+  !yCL(V1+f [3^  ’ sy^5 


which  can  be  solved  using  direct  methods. 


-u 


11) 

ip  3y ' 


( A- 1 4 ) 
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APPENDIX  B.  THE  REGRESSION  EQUATIONS 


A.  FORWARD  INTEGRATION  (1975-78  DATA) 

The  regression  equations  which  follow  were  produced  by 
SPSS  using  only  predictors  derived  from  forward  integration 
of  the  TCM.  The  data  sample  consisted  of  90  cases,  or  runs 
of  the  TCM. 


1 . TCM  Integrated  to  72  HR 

% 

VARIANCE 

REGRESSION  EQUATIONS  EXPLAINED 

DXFR12=-110.7609+  1 . 5855 ( DX1224 )+  0 . 91  54 ( 0 Y2436  ) 52.4% 

- 13. 91 66 (VX0024)+  4 . 787 5 ( X XLAT ) 

- 0. 3990(DY0012)-  0 . 231  9 ( JULDA Y ) 

+ 2 . 1836 ( VX6372  ) 

DYER12*  33.3375+18. 1010 ( V Y1224 )-  0 . 56 74 ( 0Y0036 ) 60.8% 

+ 4 . 4547 ( XXLAT ) - 1.1054(0X3648) 

+ 17.0027(VX3660)-  3 . 41 88( VXS072 ) 

- 0.3610(DY6072)-  0 . 2 0 33 ( JULDAY ) 

DXER24--261.2548+  0.0542(0X2435  )+  2 . 220 7 ( DY2436  ) 51  . 5% 

+ 1 1 . 22 15 ( XXLAT ) - 0 . 2442 ( DY4860 ) 

+ 13.0592(VY1224)-19.3594(VX0012) 

- 20.6092(VY00  36)-  0 . 4Q07( JULDAY  ) 

+ 245. 7376(VX 0036)+  6 . 1 574 ( VX60  72  ) 
-289.1865(VX0048)+67.2472(VX3648) 

DYER24=  371.6372+  4 . 70 59 ( XXLAT ) - 63 . 0 716 ( VX6072  ) 54.9% 

+ 6.3397(VY4872)-10.9032(VY0012) 

- 0.0693(DX1224)+18.3609(VX2436) 

- 1 . 4547 ( XXLON ) - 4.9576(0X0048) 

+326.2109(VX0072)+  1 . 51 70 ( DY0024 ) 

- 41.8338(VY0072)-49.6258(VX4860) 

DXER36=-344.7210+  2 . 2 548 ( DY2436 )+  0.7540(0X1224)  45.0% 

+ 15.356  5(XXLAT)  +20 . 9 332 ( V Y 1224  ) 

+ 0.8946(0X6072)-  1.0326(0X3643) 

- 0. 5088( JULDAY)-27.7637(VY0036) 

- 24.1194(VX0012)+27.4070(VX0036) 

DYER36*  260.5690-  0.9549(0X0024)-  2 . 92 55 ( DY2436 ) 55.2% 

- 1.3287(0X3648)-  0 . 5772 ( D Y60 72 ) 

- 15.52  79(VX6072)+  0 . 5606 ( DX0072  ) 

+ 3. 5077 ( DY0036 ) 
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43.8% 


DXER43*-  37.0420-  3 . 6 91 5 ( V Y2448)+  2.9949(0X1236) 

- 96.8508(VX0048)+16.3630(XXLAT) 

+ 1 . 8046 ( DY1 236 )+  8.7181 (VX5072) 

+ 1.9662(0X2436)-  2 . 4988 ( XXLON ) 

DYER483  421.5829-  6 . 0309 ( V X001 2 ) - 2 . 2291 ( DY2436 ) 

- 20.5254(VX00  72)-  0 . 50 54 ( 0Y6072) 

DXER5D*  97.7600-  9 . 8039 ( V Y2448) +16 . 2227( VX6072 ) 
+ 21 . 3398( XXLAT ) + 6 . 1 324 ( DX2436 ) 

- 59 . 0393 ( VX2448)+  1 . 9 7 5 3 ( D Y 1 2 36 ) 

- 22 . 6283 ( VXOO  12  )-  3 . 7782 ( XXLON ) 

+ 0 . 7277( DY6072) 

OYERSO-  426. 2759-14. 0207(  VX0012)-'  3..  9702  ( DY2436 ) 

- 0.1327(0X2448)-16.1361(VY4872) 

+ 0 . 9047 ( DY0036 )-  8 . 6688 ( VX60 72  ) 

DXER72*  332.5987+  7 . 1491  ( D Y24 36  ) + 20 . 2 1 35 ( V X6072  ) 
+ 33 . 3952 ( XXLAT ) - 5 . 16 16 ( XXLON  ) 

- 9.6322(VY0312)-  1 . 0394 ( JULDA Y ) 

- 31 . 3262 ( VY4860 )+  3 . 7294 ( 0X2436  ) 

- 41 .6253( VX0048) 

OYER 72*  513.8513-19. 1324(VX0012)-34.7630(VY4872) 

- 1 1 . 1 924 ( VX60 72 ) - 3 . 0878 ( DY2436 ) 

+ 32.4353(VY0024) 

2 . TCM  Integrated  to  60  HR 

REGRESSION  EQUATIONS 

DXER12=-197.0123+  1.6448(0X1224)+  1 . 1266 ( DY24 36  ) 

- 11 . 9798(VX0024)+  4 . 3382 ( XXLAT ) 

- 0.3577(DY0012)-  0 . 21 52 ( JULDA Y ) 

- 0 . 2897 ( DY4360 ) 

DYER12*  18. 0425  + 15. 8288( VY1224  )-  0 . 5743 ( DY0036  ) 
+ 3. 5755 ( XXLAT)  - 0 . 9042 ( DX3648) 

+ 8 . 4002 ( VX3560  ) 

DXER24  = -294.5950+  0.5243(0X2436)+  2 . 3803 ( 0Y2436  ) 
+ 10.6428(XXLAT)  - 0 . 6 9 73 ( DY4860 ) 

+ 15. 3517 (VY1224)— 18.4575(VX0012) 

- 16 . 437 3 ( V YOO 36  ) - 0 . 3287 ( J'JLDAY  ) 

+ 24 . 486  5 ( VX0036 ) 

DYER24*  315.0132+  3. 8023( XXLAT)  - 1 . 2271 ( DX3648) 

- 1 . 0044 ( 0Y2436 ) + 7 . 6286 ( VX3660  ) 

- 0 . 9660 ( XXLON  ) 


55.5% 

37.1% 

61.5% 

33.3% 

67.5% 

% 

VARIANCE 

EXPLAINED 

52.0% 

52.8% 

47.7% 

37.9% 
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41.1% 


0XER35--401 .1934+  2..  346 1 ( 0Y2436  ) + 2 . 7807 ( DXl 224  ) 
+ 14.181 2(XXLAT)  +25. 5785 ( VY1224) 

- 29 .8364 ( V Y0060 ) + 2.0807(0X2436) 

- 70 . 5401 ( VX0048)+  0.9644(0X2448) 

- 0 . 3235 ( JULDAY  ) 


DYER36=  323.7302-  0 . 301 0 ( DX0024  ) - 2 . 2997 ( DY2436  ) 46.5% 

- 0.9578( DX3648) 

DYER43=  11. 4836-10. 2403(VY2448)+  2 . 9591 ( DX1236 ) 42.3% 

-127.1613(VX0048)+15.9606(XXLAT) 

+ 2 . 0626 ( DY1236 )+  0.7195(0X0060) 

+ 1.8145(0X2436)-  2 . 7430 ( XXLON  ) 

DYER48=  417.1723-  6 . 5997 ( VX0012 ) - 2 . 4456 ( DY2436  ) 52.1% 

- 20.9459(VX0060) 


DXER60=  121. 4777  + 19. 3487 ( VY2448 ) +6 5 . 1 299 ( VX2436  ) 33.3% 

+ 20 .7381 ( XXLAT ) + 2 . 2910 ( DY1 236  ) 

- 23.5428(VX0012)-30.3043(VX2448) 

- 4.0537(XXL0N)  -24 . 8873 ( VY 3660  ) 

DYER50=  456. 1107-12. 666 9 ( VX001 2 ) - 3 . 8492 ( DY2436 ) 59.1% 

- 0 . 6316 ( DX2448)+23. 5083( VY0036 ) 

- 0 . 7763 ( 0 Y4860 ) 


3.  TCM  Integrated  to  48  HR 


REGRESSION  EQUATIONS 


VARIANCE 

EXPLAINED 


DXER123 -104.2041+  1 . 7146 (0X1224 )+. 0 .8720 ( DY2436 ) 50.7% 

- 14.3293(VX0024)+  3 . 648 7 ( XXLAT) 

- 0.3320(DY0012)-  0 . 1 56 2 ( JULDA Y ) 


DYER  1 2=  22. 7200+18. 7975 ( VY 1 224 ) - 0 . 5475 ( DY0036 ) 59.6% 

+ 3 . 8 7 7 3 ( XXLAT ) - 0.0330(  0X3648) 

+ 15.9946(VX2436)-21.5329(VX0048) 

- 0.3737(0Y3648) 


DXER24*-313.4117+  1.5624(0X2436)+  1 . 0949 ( DY2436 ) 46.6% 

+ 7 . 22 52 ( XXLAT)  +1 1 . 46 32 ( DX 1224  ) 

-518. 1832 (VX0048)+  0 . 7746 ( DY1 224 ) 

+ 10.5539(0X2448)+  9.8502(0X0012) 


DYER243  332.6320+  4 . 776 1 ( XXLAT ) - 0.8977(0X3648)  45.4% 

- 1 . 6 71 5 ( 0 Y2436 ) - 0 . 5969 ( DY0048) 

+ 29.4431 (VY1236)-  1 . 056 1 ( XXLON  ) 

- 6.7302(VX0012)+15.0708(VX2436) 

+ 0.4376(DY0012)-  8 . 6 575 ( VX1 236  ) 
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0XER36S-449.5619+  1 . 2046 ( DY2436 ) +14 . 0084 ( DX1224 ) 40.0% 

+ 1 1 . 1 36  9 ( XXLAT ) + 1 6 . 81 31 ( V Y 1 224 ) 

+ 1.9548(0X24 36 )-513.9317(VX0048) 

+ 12.3598(DX2448)+  11.3264(0X0012) 

DYER35®  323.7302-  0 . 301 0 ( DX0024 ) - 2 . 2997 ( DY2436  ) 46.6% 

- 0.9573(0X3548) 

DXER48=-512. 1396-13. 7846 ( VY 2448 ) + 3 . 2 353 ( DX1 236 ) 39.8% 

- 91 . 21 84 ( VX0048 ) + 1 3 . 7796 ( XXLAT ) 

+ 1 . 84 9 3 ( DY1236 ) + 1 . 721 7 ( DX24 36  ) 

DYER48=  415.0658-  6 . 31 1 2 ( VX001 2 ) - 1 . 86 54 ( DY2436 ) 52.0% 

- 0.7522(0X2448)-  0 . 4474( DY2448) 

% 

VARIANCE 
EXPLAINED 

DXER12=-104.2041+  1 . 71 46 ( DX1 224 ) + 0 . 8720 ( DY2436 ) 50.7% 

- 14. 3293(VX0024)+  3 . 6 487 ( XXLAT ) 

- 0. 3320(DY0012)-  0 . 1 562 ( JULDAY  ) 

DYER1 2 = - 5. 4141  + 18. 6 1 96 ( V Y 1 224 ) - 0 . 5781  ( DY0036  ) 51.9% 

+ 4 . 4260 ( XXLAT ) - 0.2211(0X1224) 

+ 12.1161 (VX2436)-  0.3848(0X0036) 

DXER24--290.9098+  0 . 9820 ( DX24 36 ) + 1 . 8422 ( DY2436 ) 44.2% 

+ 7 . 4308 ( XXLAT ) + 1.3908(0X1224) 

- 20. 7300(VX0024)+  1 2 . 9281 ( V Y 1 224 ) 

- 14.9940(VY0036) 

DYER24=  331.9491+  5 . 4103 ( XXLAT ) - 0 . 5106 ( DY2436 ) 31.4% 

- 0.6443(DX1224)-  4 . 7078( VY0012  ) 

- 1 . 1668(XXL0N) 

DXER36=-431 .3265+  2 . 2294 ( DY2436  ) + 2 . 2 999 ( DX1 224  ) 37.8% 

+ 11.7478(XXLAT)  + 2 1 . 66 34 ( V Y 1 224 ) 

+ 2 . 4042 ( DX2436  ) - 1 . 3375 ( DX0036  ) 

- 1 9 . 5 337 ( V Y0036  ) 

DYER36=  445.8439-  0.5419(0X0024)-  2 . 7621 ( DY2436 ) 42.9% 

+ 5.6707(XXLAT)  + 1 9 . 5 16 8 ( V Y 1 236 ) 

- 1 .4233(XXL0N) 


4.  TCM  Integrated  to  36  HR 


REGRESSION  EQUATIONS 
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The  regression  equations  which  follow  were  produced  by 
SPSS  using  predictors  derived  from  forward  and  backward  in- 
tegration of  the  TCM.  The  data  sample  consisted  of  33  cases, 
or  runs  of  theTCM  in  1977  and  1978.  Predictors  in  the  equa- 
tions are  listed  in  order  of  decreasing  significance. 

1 . TCM  Integrated  to  72  HR,  and  -36  HR  % 

' “ VARIANCE 

REGRESSION  EQUATIONS  ' EXPLAINED 

DXER 1 2=  133.7858-  1 . 1 4 34 ( XXLAT ) +49 . 0 16 2 ( VX 1 2 36  ) 85. 1% 

- 0 . 8544(  JULDAY  ) - 1 . 3384  ( DX0036  ) 

+ 0 . 8434 ( DYM200 ) + 0 . 2349 ( DXM6M2  ) 

+ 8 . 1 974 ( V YM5M4 ) - 1 . 95 32 ( VX6 0 72 ) 

+ 0.2739( DY5072) 

DYER12=-323. 6042-16. 7634(VY0012)+11 . 781 2 (XXLAT)  86.5 % 

+ 38.0056(VY0060)-  2 . 31 18( VYOO 72  ) 

- 0 . 7205 ( DY2448 ) + 21 . 8282 ( VX2436  ) 

- 34 . 35 95  ( VXOO 36  ) + 0 . 6 22 3 ( JUL DA Y ) 

- 0 . 35 79 ( DXM6M4 ) + 0 . 2499 ( 0XM4M2  ) 

DXER24=1499.882+  4 . 0 720 ( V YM200 ) - 1 . 01 76 ( DYM6M4)  88.9% 

- 2 . 7 329 ( JULDA  Y ) + 8 . 9808( VX1236  ) 

+ 0 . 0 302 ( DXM6M2 ) - 4 . 7254 ( X XLON  ) 

+ 1 . 35 1 3 ( DY2436 ) + 19 . 7881 ( VXM400 ) 

DYER24=  132. 9015  + 13. 1693(XXLAT)  - 0 . 4689 ( DY0024  ) 77.3% 

- 12 .0194(VX6072)-  1 . 2499 ( XXLON ) 

+ 10.8862(VX2436)-  0 . 31 6 1 ( DXM6M4 ) 

+ 1 2 . 57  33  ( VY1 224 ) - 0 . 5091  ( DY3648) 

+ 0 . 5794 ( DX4872) - 1 2 . 1 879 ( V X00 36  ) 

DXER36  = 2313. 2390  + 14. 1 02 5 ( V YM200  ) - 1 . 46 77 ( DYM6M4  ) 89.6% 

- 4 . 9703 ( JULDAY ) + 1 3 . 8301  ( VX1236  ) 

+ 1 8 . 86 0 7 ( VXM4M2 ) - 4 . 884 1 ( XXLON  ) 

- 12 . 2852 ( XXLAT) 

DYER36=  208. 2644  + 13. 0715(XXLAT)  - 2 . 1 86 3 ( X XLON  ) 74.9% 

- 2 . 3350 ( VX4872 ) + 0 . 2046 ( VY2436 ) 

- 0 . 7327 ( DXM6M4 ) + 1 2 . 36 56 ( VYM6M4 ) 

+ 0.9038( JULDAY) 
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DXER48=  3766 
+ 2 
- 10 
+ 24 
+ 27 


9010+  3 . 96  98 ( V YM200  ) - 
92  78 ( VX48  72 ) - 27.2623 

1 248 ( VY0048 ) - 18.3859 

7539 ( VXM4M2  )-  7.2417 

7765( VX1236 )-  1.3987 


7 . 5053 ( JULDAY) 
( XXLAT ) 
(VYM6M4) 

(XXLON ) 

( DX3648) 


92.4% 


DYER43=  587. 

- 5. 

- 4. 
+ 8. 
+ 2. 
+ 20. 


8164-  1 . 24 70 ( DXM6M2  ) - 1 
9702(VX4872)+  91.4204 
3418( XXLON)  + 0.2002 
8435 ( XXLAT ) + 1.8326 
0270  ( DXM4M2)  - 1.424.5 
3420 ( V YM6M4 ) - 1.3766 


3 . 4294 ( VY2436  ) 

(VYM424) 

(0X3648) 

(JULDAY) 

(DYM212) 

(-  0 Y M 2 0 0 ) 


81.0% 


DXER60  = 4803 .057  +1 9 . 06 39 ( V YM200 ) - 10 . 5486 ( JULDA Y ) 88.6% 

+ 8.5565(VX4872)-  42 . 8314 ( XXLAT ) 

- 46 . 1 492 ( VY0024 ) - 6 . 9290 ( XXLON ) 

+ 2 . 3348 ( DXM4M2 ) - 1 . 66 72 ( 0 YM6M4 ) 


DYER60=  364.9735+  0 . 6 71 1 ( DXM600 ) + 4 . 1587 ( JULDAY  ) 80.9% 

+ 21 .0652( XXLAT)  + 7 5 . 96 30 ( V Y0024 ) 

+ 2.0039(DYM6M2)-  7 . 1 881 ( XXLON ) 

- 2 . 7226 ( DXM6M4  ) - 2 . 060 2 ( DYM2 1 2 ) 


DXER72  = 3715 .035  + 87 . 8945 ( VYM200  ) - 6 . 6852 ( DYM6M4 ) 95.7% 

- 12. 7603( JULDAY)-  115 . 0140( VY0048) 

+ 42 . 7842  ( V XM6M2 ) - 39 . 1 6 36 ( XXLAT ) 

- 0.8422(0X6072) 


DYER72*1214. 7320-26. 6 951 ( V YM200  ) - 7 . 6 1 89 ( XXLON  ) 93.2% 

- 2 . 7921 ( DXM6M4)+  6 2 . 87 56 ( V YM6M4 ) 

+ 38.9121(VY0024)+  3. 1 0 78 ( JULDAY  ) 

+ 1 .8961 (DY3648)+  2 1 . 5 76 3 ( V X 3648  ) 

- 23 . 6 336 ( VX001 2 ) 


2. 


TCM  Integrated  to  60  HR,  and  -36  HR 


REGRESSION  EQUATIONS 


% 

VARIANCE 

EXPLAINED 


DXER1 2=  167.5636-  1 . 1 6 79 ( XXLAT)  +42 . 0954 ( VX1 2 36  ) 

32.2% 

- 0.9556( JULDAY)- 

1.1599(0X0036) 

+ 0 . 8786 ( DYM200  ) + 

- 5 . 6 377 ( V YM6M4 ) 

0 . 2185 ( 0XM6M2 ) 

DYER12=-344. 2635-14. 9407( 

VY0012)+12.0804(XXLAT) 

86.4% 

+ 26 . 545  5 ( VY0060  ) + 

1.6737(0X2436) 

- 26 . 5420 ( VXOO 36  ) - 

0.6601 (DY3648 ) 

+ 0 .6473( JULOAY  )- 

+ 4 . 9116 ( VXM200) 

0 . 300 1 ( DXM6M4 ) 
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DXER24* 1499. 8820+  4 . 0 720 ( V YM200 ) - 1 . 01 76 ( DYM6M4 ) 

- 2 . 7 3 2 9 ( JULDAY )+  8.9808(VX1236) 

+ 0 . 0302 ( DXM6M2 ) - 4 . 7254 ( XXLON ) 

+ 1 . 35 1 3 ( 0Y2436 ) + 19 . 7881 ( VXM400  ) 

i 

88.9% 

1 

DYER24  = - 56. 4323+12. 9000(XXLAT)  - 0 . 3905 ( 0Y002 1) 
+ 9 . 0 748 ( VY4860 ) - 0 . 3743 ( 0XM6H4 ) 

<1  + 0.5430( JULDAY)-  C . 96 95 ( XXLON ) 

72.5% 

; ' 

j DXER36  = 2312. 2390  + 14. 1025( VYM200) - 1 . 4677 ( DYM6M4) 

- 4. 9 7 0 3 ( JULDAY  )+  1 3 . 8301  ( VX1 236  ) 

+ 1 8. 860  7 ( VXM4M2 ) - 4 . 8841 ( XXLON ) 

- 1 2 . 2852 ( XXLAT ) 

89.6% 

DYER36=  518. 1446  + 14. 5509(XXLAT)  - 3'.  906  7 ( XXLON  ) 

- 4 . 4 0 7 7 ( VX4860 )+  1 . 0088 ( JULDAY ) 

+ 1 . 25  97 ( DYM6M2 ) - 1 . 6445 ( DXM6M4 ) 

+ 19.1652 ( VXM600 ) + 9 . 26 26 ( V Y 1 224 ) 

| - 1 3 . 881 3 ( V YM400 ) 

78.4% 

DXER48-3901.1060+  0 . 1432 ( VYM200  ) - 7 . 1 76 5 ( JULDA Y ) 
j - 26 . 0746 ( XXLAT ) + 33. 2885 ( VX1236 ) 

- 8. 5896 ( XXLON)  + 0 . 4893( D Y243o ) 

If  + 2 . 0444  ( DXM4M2  ) - 1 . 36  76  ( DX  3648 ) 

J * - 1 . 3292 ( DYM6M4 ) 

92.2% 

DYER48=1228.8470-  3 . 0296 ( DXM6M2 ) - 9 . 980 1 ( V Y24 36  ) 

+ 56 . 0325 ( VYM6M2 ) + 1 .42 10 ( JULDAY ) 

+ 50 . 1 396  ( VY0024 ) - 4 7 . 0803 ( V YM21 2 ) 

- 5 . 8375 ( XXLON ) + 2 . 551 8( DXM4M2 ) 

+ 43 . 1 560 ( VXM600 ) - 7 . 6 1 43 ( VX2436 ) 

81.4% 

DXER60  = 5564. 0920-  7 . 204 1 ( VYM200 ) - 10 . 5 756 ( JU LDA Y ) 

- 55 ,6857( XXLAT ) - 1 247 . 4050 ( VY 0024 ) 

1 + 21 . 4406 ( VX4860 ) - 9 . 8929 ( XXLON ) 

+ 1833.842  ( VYOO  36  ) + 2 . 3 39 2 ( DXM4M 2 ) 

-620.6515(VY2436)+  26 . 8308 ( VY2448 ) 

90.5% 

DYER603  364.9736+  0 . 6 71 1 ( DXM600 ) + 4 . 1 587 ( JULDA Y ) 

+ 21 . 06  52 ( XXLAT ) + 75 . 96 30 ( VY0024 ) 

+ 2 . 0889 ( D YM6M2 ) - 7 . 1 881 ( XXLON ) 

- 2 . 7226  ( DXM6M4  ) - 2 . 0602 ( DYM2 12  ) 

80.9% 

I 84 

? 

Timr- 

u 

3.  TCM  Integrated  to  48  NR,  and  -36  HR 


REGRESSION  EQUATIONS 


DXER12*  167.5636-  1 . 1 6 7 8 ( XXLAT)  +42 . 0954 ( VX1236 ) 

- 0 . 9556 ( JULOAY ) - 1 . 1 5 99 ( DX0036 ) 

+ 0 . 8736 ( DYM200  ) + 0 . 2184 ( 0XM6M2  ) 

- 5 . 6 377 ( V YM6M4 ) 

DYER12=-283.4232-  8 . 4050 ( V Y001 2 ) + 9 . 2910 ( XXLAT) 

+ 12.9974(VY1236)+  1 . 2962 ( DX2436  ) 

- 17.6267(VX0036)+  0 . 6 355 ( JULDA Y ) 

+ 5 . 6809 ( VYM6M4) - 0 . 2878 ( DXM6M4 ) 

DXER24*1499.8820+  4 . 0 720 ( V YM200 ) - ' 1 . 01 76 ( 0YM6M4  ) 

- 2.732  9 ( JULOA Y ) + 8 . 9808 ( VX 12 36  ) 

+ 0 . 0 302 ( DXM6M2 ) - 4 . 7254 ( XXLON ) 

+ 1 ,3513( DY2436)+19. 7381 (VXM400) 

DYER24=  129.5675+11. 7946(XXLAT)  - 0 . 4375 ( DY0024 ) 
+ 1 0 . 8043 ( V Y 1 236 ) - 0. 5345 ( DXM6M4) 

+ 7 . 94 36  ( V YM6M4 ) - 1 . 4 770 ( XXLON  ) 

+ 0.5032( JULOAY) 

DXER36=2313. 2390+14. 1025(VYM200)-  1 . 46 77 ( DYM6M4 ) 

- 4.9703(JULDAY ) + 13.8301  (VX1236) 

+ 1 8 . 8607 ( VXM4M2 ) - 4 . 8841 ( XXLON ) 

- 12.2852(XXLAT) 

D YER36=  130. 4497+13. 4500(XXLAT)  - 1 . 93 19 ( XXLON ) 

- 0. 8659(0 XM6M4)+13.7236(VYM6M4) 

+ 0.9130 (JULOAY ) 

DXER48-3901.1060+  0 . 1 4 32 ( VYM200  ) - 7 . 1 76 5 ( JULDA Y ) 

- 26 . 0746 ( XXLAT ) +33 . 2885 ( VX1 2 36 ) 

- 8.5896  (XXL0N)  + 0 . 4893 ( D Y2436  ) 

+ 2 . 0444 ( DXM4M2  ) - 1.36  76(  0X3648) 

- 1 . 3292 ( DYM6M4 ) 


DYER48=1228.847  - 3 . 0296 ( DXM6M2  ) - 9.9801 
+ 56 . 0 325 ( V YM6M2 ) + 1 . 42 1 0 ( JULDA Y ) 
+ 50.1396(VY0024)-47.0803(VYM212) 
- 5 . 8375 ( XXLON ) + 2 . 5 5 1 8 ( DXM4M2 ) 

+ 43 . 1 560 ( VXM600 ) - 7 . 6 143 ( VX2436  ) 


VY2436) 


% 

VARIANCE 

EXPLAINED 

82.2% 

85.5% 

88.9% 

75.9% 

89.6% 

74.5% 

92.2% 

81.4% 


DXER12=  167.5686-  1 . 1679 ( XXLAT)  +42 . 0954 ( VX1 236 ) 

- 0.9556(JULDAY)-  1.1599(0X0036) 

+ 0 . 8786 ( DYM200 ) + 0 . 2185 ( DXM6M2) 

- 5 . 6 377 ( VYM5M4 ) 

DYER12=-283.4232-  8 . 4049 ( V Y00 1 2 ) + 9 . 2910 ( XXLAT) 

+ 12 . 9974 ( V Y1 236 ) + 1.2962(0X2436) 

- 17.6267(VX0036)+  0 . 6 355 ( JULDA Y ) 

+ 5.6809 ( V YM6M4) - 0 . 2 878 ( DXH6M4) 

DXER24=1499.8820+  4 . 0720 ( VYM200 ) - 1 . 0 1 76 ( DYM6M4  ) 

- 2 . 7329 ( JUL DA  Y ) + 8 . 9808 ( VX 1 236 ) 

+ 0 . 0 3 0 2 ( DXM6M2 ) - 4 . 7254 ( XXLON ) 

+ 1 . 3513(DY2436 )+19 . 7881 ( VXM400) 

DY ER24=  129. 5675  + 11. 7946(XXLAT)  - 0 . 4375 ( DY0024 ) 
+ 10.8043(VY1236)-  0 . 5345 ( DXM6M4 ) 

+ 7 . 94  36 ( V YMSM4 ) - 1 . 5 770 ( XXLON  ) 

+ 0.5032( JULDAY) 

DXER36=2313. 2390+14. 1025(VYM2Q0)-  1 . 46 77 ( DYM6M4 ) 

- 4.9703(JULDAY)+13.8301(VX1236) 

+ 1 8 . 8607  ( V XM4M2) - 4 . 8341 ( XXLON  ) 

- 12 . 2852 ( XXLAT ) 

D YER36=  130. 4497+13. 4500 ( XXLAT)  - 1 . 93 1 9 ( XXLON ) 

- 0.8659(DXM6M4)+13.7236(VYM6M4) 

+ 0 . 91 30 ( JUL  DA Y ) 


% 

VARIANCE 

EXPLAINED 

82.2* 


85.5% 


88.9% 


75.9% 


89.6% 


74.5% 
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